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Abstract —  In  this  paper,  we  investigated  different  methods 
for  blind  Doppler  shift  estimation  and  compensation  for  a 
single  carrier  in  underwater  acoustic  wireless  sensor  networks. 
We  analyzed  the  data  collected  from  our  experiments  using 
non-data  aided  (blind)  techniques  such  as  Power  Spectrum 
Analysis,  Autocorrelation,  and  Squaring  Time  Phase  Recovery 
methods  in  order  to  estimate  the  Doppler  shift  in  collaborative 
distributed  underwater  sensor  networks.  Detailed 
experimental  and  simulated  results  based  on  second  order 
cyclostationary  features  of  the  received  signals  are  presented. 
Keywords:  Blind  Doppler  Shift  Estimation,  Underwater 
Communication,  Autocorrelation,  Power  Spectral  Density 
(PSD),  Periodogram. 

i.  Introduction 

Doppler  shift  estimation  and  detection  for  target 
localization  and  tracking  in  underwater  wireless 
communication  (UWC)  has  been  a  major  topic  of  research 
and  investigation  due  to  the  increasing  use  of  aquatic 
channels  [2],  [4],  [12],  and  [15].  The  need  for  Doppler  shift 
estimation  in  UWC  exists  mostly  for  real  time  remote 
control  monitoring  of  oceanic  activities:  environmental 
monitoring,  scientific  data  collection,  tracking,  and  locating 
objects. 

There  are  several  primary  obstacles  for  reliable 
communication  in  underwater  environment,  including  time- 
varying  multipath,  fading,  low  sound  speed  and  noise.  The 
sound  speed  underwater  is  about  1500  m/s.  The  low  sound 
speed  and  relative  high  platform  speed  result  in  Doppler 
shifts  several  times  those  encountered  in  radio  transmission 
[14],  The  ratio  of  platform  speed  to  sound  propagation  speed 
is  large  enough  to  cause  time  compression  or  expansion  of 
the  symbol  pulse  itself  [4],  [13],  The  motion-induced  pulse 
compression  or  expansion  makes  symbol  synchronization  of 
equal  importance  to  carrier  frequency  identification,  thus 
becoming  a  major  constraint  of  mobile  UWC.  Current  UWC 
methods  are  mostly  used  for  low  Doppler  environment  [4] . 

Due  to  complexities  of  the  underwater  channel,  such  as 
multi-path  propagation,  time  variations,  small  available 
bandwidth  and  strong  signal  attenuation,  (especially  over 
long  ranges)  various  factors  of  a  communication  system 
such  as  data  rate,  symbol  synchronization,  carrier  phase 


recovery  and  the  speed  of  propagation  are  compromised  [3]. 
As  a  consequence,  at  the  digital  receiver  end,  different 
levels  of  synchronization:  carrier  recovery,  frame 
synchronization,  symbol  and  data  bit  timing  recovery  are 
strongly  affected  [9],  [10]  and  [14].  Therefore,  continuous 
time  communications  between  rapidly  moving  platforms  is 
the  driver  for  new  robust  methods  for  blind  synchronization 
(non-data  aided)  techniques  able  to  track  large  and  variable 
Doppler  shifts. 

In  order  to  estimate  and  compensate  for  the  Doppler 
shift,  different  coherent  and  blind  methods  have  been 
implemented  for  pseudorandom  sequence  estimation.  The 
most  popular  can  be  categorized  as  eigenanalysis  based  on 
subspace  methods,  and  exploiting  the  spectral  characteristics 
in  cyclostationarity  signals.  In  [2],  [4]  it  has  been  proven 
that  the  spectral  characteristics  of  cyclostationarity  signals  is 
computationally  less  complex  and  more  robust  for 
pseudorandom  (PN)  sequence  estimation  then  subspace 
methods.  Therefore,  we  have  chosen  the  experimental 
approach  to  estimate  the  Doppler  shift  using  blind  spectral 
estimation  methods  for  data  modulated  by  PN  sequences. 

In  [11],  spectral  correlation  based  signal  detection  has 
been  proposed.  The  spectral  correlation  theory  in  [11]  is 
used  to  calculate  spectral  correlation  function  and  it  could 
be  used  for  Doppler  shift  estimation.  However,  in  this 
method,  the  received  baseband  signal  is  no  longer 
orthogonal  to  the  transmitted  //(-sequences.  As  a  result,  in 
the  frequency  domain,  it  is  very  hard  to  read  the 
instantaneous  Doppler  shift  due  to  fading. 

In  [6]  the  concept  of  passive  signal  detection  is  carried 
forward  to  active  signal  detection  using  a  Dopplergram  and 
an  ambiguity  function  has  been  used  to  determine  Doppler 
shift  for  ///-sequence  modulation.  However,  in  both  [6]  and 
[11],  only  one  method  of  modulation  (///-sequence 
modulation  technique)  has  been  used.  In  addition,  in  [4],  the 
spectral  correlation  function  is  modified  to  the  spectrum 
coherence  function  to  estimate  carrier  frequency  and  symbol 
rate  estimation.  It  was  assumed  that  in  underwater 
communication,  channel  characteristics  vary  quickly  and  the 
signal  parameters  vary  quickly  as  well. 

In  this  paper,  we  conducted  experiments  for  underwater 
acoustic  wireless  communication  using  a  pair  of  SAM-1 
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sensors  provided  by  Desert  Star  Systems.  We  used  universal 
asynchronous  receiver  transmitter  protocol  for  serial  data 
communication.  MATLAB  software  was  used  to  send  data 
via  the  serial  port  for  transmission  and  acquire  data  from  a 
sensor  into  a  personal  computer  (PC).  We  analyzed  the  data 
collected  from  our  experiments  using  non-data  aided 
techniques  such  as  Power  Spectrum  analysis, 
Autocorrelation  and  Squaring  Time  Phase  Recovery  (Oerder 
&  Meyr)  [1]  methods  in  order  to  estimate  Doppler  shift  in 
collaborative  distributed  underwater  sensor  networks.  In  our 
study,  the  sensors  were  half-duplex,  and  therefore  could 
only  transmit  or  receive  at  a  given  time. 

We  improved  the  MATLAB  code  for  serial  data 
communication  for  acoustic  sensors  (SAM-1)  provided  by 
Desert  Star  Systems.  We  transmitted  original  and  modulated 
52  m-sequences  each  of  length  1023  bits,  via  sensors  in  an 
acoustic  prototype  environment  and  at  the  receiver  end  we 
analyzed  the  received  signal  using  spectral  analysis  and  the 
Oerder  and  Meyr  method.  The  size  of  the  baseband 
transmitted  covariance  matrix  is  106392x106392.  The 
received  signal  correlation  matrix  is  1x106392  for  the 
original  m-sequences. 

The  rest  of  the  paper  is  organized  as  follows:  In  Section 
II  we  describe  the  theoretical  background.  Section  III 
presents  the  experimental  model.  Section  IV  contains  our 
simulated  and  experimental  results.  Conclusion  and 
acknowledgements  are  drawn  in  section  V. 

II.  Theoretical  Background 
1)  Doppler  Effect 

Generally,  change  in  the  frequency  of  an  emitted  wave 
caused  by  the  motion  of  an  emitted  source  relative  to 
observer  or  vice  versa  is  defined  as  Doppler  shift  or  effect. 
For  a  communication  system,  the  received  signal  at  the 
receiver  end  can  be  characterized  as: 

r(t)  =  9I(Xn=-oo  dn  git-  nT )  el2nfct}  +  n(t)  (1) 

where,  91  represents  real  part  of  signal,  n(t)  is  the  AWGN 
with  zero  mean  circular  complex  white  Gaussian  process 
statistically  independent  signal,  fc  is  the  carrier  frequency, 
dn  is  the  transmitted  data  symbol  in  the  time  interval 
nT  <  t  <  fn  +  1)7,  T  is  data  symbol  duration  and,  g(t)  is 
the  convolution  product  of  impulse  response  of  pulse 
shaping  filter  g£(t),  channel  impulse  response  h(t),  and 
receiver  filter  impulse  response  gR  (t)  [12]. 

In  wide  band  cases,  due  to  the  Doppler  effect,  signal 
carrier  frequency  suffers  frequency  scaling  /c(l  —  Sf)  and 
the  received  baseband  signal  undergoes  time  scaling 
t(l  —  8d ),  so  the  received  signal  at  the  receiver  input  is 
given  in  [2]  and  [16]  as: 

r(t)  = 

9*{l£=-ood„  giti  1  -  So)  -  n7V-/c(i-^}  +  n(t) 

(2) 


where  8d  is  the  relative  Doppler  shift.  This  relation  is  valid 
for  both  wireless  communication  systems  and  UWC 
systems.  But,  the  Doppler  shift  in  underwater 
communication  is  very  high.  For  underwater  communication 
with  multipath  and  fading,  the  received  signal  can  be  given 
as: 

rit,8d )  =  ^a^-ooAf^^dngm  -  Sd)  - 
-nT]ej2n[^+A^tCl-s^}  +  n(t) 

(3) 

where  Aft)  is  the  fading  gain  of  Ith  path,  L  is  the  number 
of  multipath  component,  and  (t)  is  the  phase  offset  due  to 
the  channel  on  Ith  path. 

Different  coherent  techniques  have  been  used  to 
estimate  Doppler  shift  but  the  proposed  algorithms  used  to 
find  cyclic  frequency  offset  are  more  susceptible  to  ISI 
(Inter  Symbol  Interference).  Thus,  non-coherent  techniques 
are  preferred  to  find  cyclic  frequency  offset  and  Doppler 
compensation  [7],  [8],  [16]  and  [17]. 

2)  Spectral  Analysis 

Most  random  processes  encountered  in  nature  arise  from 
some  periodic  phenomena.  The  random  processes  generated 
from  such  periodic  phenomena  produces  data  that  are  not 
periodic  functions  of  time,  but  their  statistical  properties 
varies  with  time.  These  kinds  of  random  processes  are 
modeled  as  wide-sense  cyclostationarity  random  processes 
and  its  features  can  be  used  in  signal  detection  and 
estimation.  Therefore  random  processes  in  this  paper  are 
considered  to  be  cyclostationary. 

Cyclostationary  analysis  is  based  on  the  fact  that 
communications  signals  are  not  accurately  described  as 
stationary,  but  rather  more  appropriately  modeled  as 
cyclostationary.  While  stationary  signals  have  statistics  that 
remain  constant  in  time,  the  statistics  of  cyclostationary 
signals  vary  periodically.  These  periodicities  occur  for 
signals  of  interest  in  well-defined  manners  due  to  underlying 
periodicities  such  as  sampling,  modulating,  multiplexing, 
and  coding. 

A  process,  sayx(t),  is  said  to  be  wide  sense 
cyclostationary  if  it’s  mean  (if (A:(t)})  and  autocorrelation 
function  (Rxit,  r))  are  periodic  with  the  same  period  T : 

E{xit  +  T)}  =  E{xit)}  (4) 

Rxit  +  T,t)  =  Rxit,x)  (5) 

In  [7],  the  Fourier  series  expansion  of  this  periodic 
autocorrelation  function  converges.  As  in  [14],  (5)  can  be 
expressed  as: 

Rxit, r)  =1.^-0,  Rlir)ei2n(T>  (6) 

n 

where  -  =  a  is  called  the  cyclic  frequency;  Rxi r)  is  the 
cyclic  autocorrelation  function  at  cyclic  frequency  a  and  is 
given  as  follows: 


Rx(j)  A  lim7._0O \\T_jj2  Rx(t.  *)e  ,2nat  dt  (7) 

Now  for  random  process,  the  Fourier  transform  of  the 
autocorrelation  will  give  the  Power  Spectral  Density  (PSD) 
and  (7)  can  be  expressed  in  so  called  cyclic  power  spectral 
density  as: 

W)  =  ZJSL.  Rax(T)e-’2*r*  (8) 

The  cyclic  PSD.  5“(/).  contains  spectral  discrete 
components  which  are  useful  to  estimate  the  Doppler  shift 
for  different  modulated  signals. 

Furthermore,  spectral  components.  S“(/0,  of  signal 
x(t)  are  the  measurement  of  cyclic  pow'er  spectral  density 
which  can  be  more  elaborated  by  the  normalized  correlation 
between  two  spectral  components  of  x(t)  at  /  +  jand  f  — 

j  frequencies  over  an  interval  of  At.  Then  the  ideal 
measurements  can  be  mathematically  expressed  as  in  [14]: 

s?(f)  = 

limr_(00  limAt_,00  ^  flt^xT(t,  f  +  |)xf  (t,  /  -  f)  dt 
2 

(9) 

where  the  finite  time  Fourier  transform  of  x(t)  over  time 
interval 

,  T  .  T  . 

t  +  -andt - is: 

2  2 

XT(t,0  =  ft-T/2  x (M)ei2nTudu  (10) 

III.  Experimental  Model 

We  conducted  experiments  in  an  indoor  water  tank,  and 
hi  an  indoor  swimming  pool.  The  tank  was  7x5x2. 5  ft3  in 
volume.  The  data  were  transferred  using  a  pan  of  acoustic 
modems  and  wrere  processed  in  MATLAB.  Also  we  used  a 
hydrophone  to  measure  sound  underwater  and  an  acoustic 
speaker  to  generate  noise.  Fig.  1  illustrates  our  prototype 
environment;  where  the  sensors  are  connected  to  the 
computer  via  a  serial  port  and  are  floating  in  the  water. 

hi  the  swimming  pool,  distance  between  two 
communicating  sensors  was  15  yards  and  the  depth  of  the 
swimming  pool  wras  lift  and  8  niches.  The  surface 
temperature  of  the  swimming  pool  was  70°  Fahrenheit  and 
the  water  was  chlorinated. 


TV.  Experimental  Results 

Our  experiment  is  based  on  serial  data  transfer.  The  data 
was  transmitted  at  4800  baud  and  8  data  bits.  We  measured 
the  noise  during  our  experiment  for  the  length  of  3.65 
minutes  sampled  at  the  rate  of  .01s  and  we  plotted  its  PSD. 

From  Fig.  2  one  can  notice  that  the  power  of  noise  is 
high  for  frequency  range  of  0  to  5  kHz.  As  the  fr  equency 


increases  we  observe  that  the  power  of  noise  is  less  than 
—11.11.86  dB.  If  wre  are  able  to  transmit  the  signal  in  the 
frequency  range  where  the  noise  power  is  relatively  less,  we 
can  obtain  better  performance.  Therefore,  we  have  chosen 
17  kHz  of  carrier  frequency  for  our  experiments.  We  could 
use  any  frequencies  above  5  kHz  and  below  40  kHz.  but  to 
be  consistent  with  results  presented  hi  [2],  [4],  we  preferred 
17  kHz. 
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Figure  1  Experimental  set  up  and  its  components 


Power  Spectral  Density 


Figure  2  PSD  of  measured  noise  m  swimming  pool 

In  our  experiments  the  transmitted  signal  had  a  carrier 
frequency  of  17  kHz  and  a  sampling  rate  of  8000 
samples/sec.  The  transmitted  spectrum  in  both  the  cases 
(tank  and  pool)  had  almost  a  flat  spectrum.  For  the  1st  case 
the  transmitted  signal  is  a  PN  sequence  of  length  1023  bits 
modulated  by  a  earner  of  17  kHz.  We  plotted  the 
autocorrelation  and  partial  PSD  of  the  received  signal  and 
observed  the  different  spectral  components  to  estimate 
Doppler  shift  for  both  tank  and  pool. 

It  is  clear  form  Fig.  3  that  based  on  the  autocorrelation 
function  of  the  received  signal  we  were  able  to  distinguish 


the  transmitted  52  ///-sequences.  However,  the  received 
signal  is  not  free  of  noise  and  some  multipath  components. 
The  dark  lines  on  Fig  .3  refer  to  the  cross-spectral 
correlation  effect.  The  bold  dense  line  on  the  autocorrelation 
of  the  transmitted  signal  is  due  to  modulation  effect.  We 
plotted  the  partial  autocorrelation  of  the  received  signal  and 
find  out  the  transmitted  ///-sequences  associated  with  it  as  in 
Fig.  4. 


Autocorrelation  Frequency  Modulated  Received  Signal 


shift.  The  Doppler  shift  was  50  Hz  at  fc  =  17kHz.  The 
partial  PSD  spectrum  in  Figure  6  shows  that  the  Doppler 
shift  is  20  Hz  at  2 fc.  Taking  average  we  conclude  that  the 
estimated  Doppler  shift  is  approximately  35  Hz. 
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Figure  5  Partial  PSD  of  a  received  signal  around  twice  the  carrier  frequency 
(results  for  water  tank) 


Figure  3  Autocorrelation  function  of  the  received  52  m  -sequences 


In  Figure  4.  we  noticed  that  the  time  length  difference 
between  each  ///-sequence  peak  is  not  1023  bits  (in  the 
transmitted  signal  it  was  exactly  1023  bits).  Form  this 
observation  we  concluded  that  the  multipath  effect  must 
have  affected  the  received  signal  during  the  period  of  1023 
milliseconds. 


Autocorrelation  ot  Frequency  Modulated  Received  Signal 
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Figure  4  Partial  autocorrelation  of  received  signals  showing  48th  to  52nd  m- 
sequences 

hi  the  case  of  the  water  tank  experiment,  we  were  able  to 
receive  same  data  that  were  sent  with  a  bit  error  rate  of 
0.0078.  Looking  at  the  partial  PSD  around  the  carrier 
frequency  and  its  period,  in  Fig.  5,  we  find  the  Doppler 
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Figure  6  PSD  of  received  signal  (results  for  swimming  pool) 


After  measuring  the  Doppler  shift  in  the  water  tank,  we 
performed  experiments  in  the  swimming  pool,  hi  the 
swimming  pool,  the  transmission  range  was  longer  (25 
yards)  and  the  environment  is  echoic,  so  we  experienced 
more  of  multipath  effect.  We  searched  for  the  local  maxima 
around  the  carrier  and  found  the  Doppler  shift  to  be  130  Hz 
at  fc  =  17kliz  and  170Hz  at  2 fc  in  average  150Hz.  Then, 
we  transmitted  the  Linear  Frequency  Modulated  ///- 
sequences  signal  and  analyzed  the  signal  using  Oerder  & 
Meyr  squaring  recovery  circuit.  Figures  7-8  show  the  partial 
PSD.  and  phase  shift. 

From  Fig.  8,  we  see  that  average  Doppler  shift  is  0.85 
and  the  instantaneous  shift  of  range  [0-2.86],  From  the 
experiments,  we  claim  that  Doppler  shift  could  be 


minimized  by  using  double  modulation  technique  instead  of 
using  either  the  LFM  or  m-sequence  technique  separately. 
This  can  be  validated  by  the  minimal  shift  obtained  in 
Figure  8. 
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Figure  7  Partial  PSD  of  LFM  52-  m-  sequences  at  2nd  period  (34  kHz) 


It  can  be  observed  that,  for  the  proposed  double 
modulation  technique,  phase  shift  is  less  than  the  one 
provided  in  [2],  The  results  provided  in  [2]  were  for  m- 
sequence  modulation  and  the  average  phase  shift  was  1.12 
for  the  same  parameters  (k=T00  symbols).  From  Fig.  8,  we 
see  that  average  Doppler  shift  is  0.85  and  the  instantaneous 
shift  of  range  [0-2.86]. 


M-l.m-sequences 


Figure  8  Estimation  of  Instantaneous  and  average  Doppler  shift  using  the 
Oerder  &  Meyr  Algorithm 


In  Fig.  9  we  focused  only  on  die  instantaneous  estimated 
Doppler  shift.  The  sampling  frequency  is  20Kz  and  the 
number  of  samples  per  symbols  is  4:  please  see  [17]  and 
[18].  We  converted  the  received  baseband  signal  in  a 
number  of  total  212784  samples  and  we  used  the  same 
number  of  K=100  symbols  per  frame.  Because  we  are  using 
a  blind  method,  we  cannot  compensate  the  Doppler  shift  for 
the  first  ///-sequence  and  w'e  obtained  a  “coarse”  estimation 
of  this. 


We  varied  the  numbers  of  symbols  per  frame  between 
50  and  100  symbols,  and  the  results  are  not  significantly 
changed.  One  possible  explanation  for  this  is  that  the 
number  of  samples  per  symbols  is  the  same  (=4),  and  from 
information  theoretic  point  of  view,  the  information 
extracted  fr  om  each  ///-sequence  is  the  same. 

It  is  an  interesting  question  to  ask.  what  could  be  the 
optimum  number  of  symbol  per  frame  and  the  minimum 
number  of  samples  per  symbol,  in  order  to  get  the  optimum 
of  performances  or  to  get  the  maximum  of  information.  This 
is  left  for  future  research  and  the  results  obtained  in  Fig.  8 
and  Fig.  9  are  a  good  start. 


Baseband  received  x(1 512784),  fs=1/20000,  4samples  x1023 


Figure  9  Instantaneous  estimation  of  Doppler  shift 


We  searched  for  a  Doppler  shift  between  [-30.  30]  Hz 
and  wre  wanted  to  plot  all  52-ambiguity  functions  for  each 
///-sequence  as  in  Fig.  11.  Please  note  a  scale  change  for 
Doppler  shift  measurement  in  Fig.  10.  In  searching  for  the 
local  maxima  or  minima  of  the  ambiguity  functions,  we 
scaled  the  domain  of  definition  for  these  functions. 
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Fig.  10  The  ambiguity  function  for  all  52  m-sequences 


hi  comparing  with  the  method  given  in  [2],  in  average, 
there  is  a  difference  of  0.5  Hz,  and  each  ///-sequence  has  a 
different  Doppler  shift.  The  explanation  for  this  is  the  fact 


that  blind  method  is  not  a  coherent  one  (data  aided)  and  it 
depends  which  sample  (out  of  four  per  symbol)  is  coming 
and  processed  first. 

The  best  comparison  is  made  directly  on  the 
Dopplergram  [5]  represented  in  Fig.ll.  From  user 
perspective  (Dopplergram)  there  is  practically  no  difference 
between  results  in  Fig.  10  and  that  those  presented  in  [5]. 
The  explanation  for  this  are  the  results  presented  in  the  Fig. 
10,  if  we  are  looking  only  at  the  maximum  of  the  ambiguity 
functions.  Based  on  Fig.10  we  can  provide  directly  the 
Dopplergram: 


Doppler  Gram 


Doppler  Shift  (Hz) 

Figure  1 1  Dopplergram  results 

During  our  experiment  wTe  had  several  issues  with  data 
and  software  compatibility.  Debugging  was  the  most 
cumbersome  part.  It  took  several  hours  to  run  the  MATLAB 
programs.  Memory  buffer  size  of  sensors  was  another  one 
challenging  issue.  We  had  problems  with  the  sensors 
because  they  transmitted  some  random  signals  by 
themselves  due  to  unknown  internal  error. 

V.  CONCLUSIONS 

In  this  paper,  we  performed  experiments  in  two  different 
types  of  environment  and  studied  the  non-data  aided 
techniques  for  Doppler  shift  estimation  in  underwater 
communication.  We  found  that  the  Doppler  shift  for  the 
large  transmission  range  wras  higher  than  for  the  lower 
transmission  range.  We  found  that  double  modulation 
technique  (combine  LFM  and  ///-sequence  for  modulation 
message  signal)  could  minimize  the  phase  shift  even  better 
than  the  modulation  technique  used  in  [3],  for  underwater 
communication. 

Based  on  our  observations,  we  conclude  that  the 
exploiting  second  order  cyclostationary  features  of  the 
received  signal  makes  it  easier  and  faster  to  estimate 
Doppler  shift  without  prior  knowledge  of  the  signal 
transmitted.  It  is  also  concluded  form  our  experiments  that 
the  proposed  methods  were  very  easy  to  implement. 
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